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Abstract 

Following the arguments presented by Mansuripur [Opt. Express 
16, 14821-14835 (2008)], we suggest a form for the macroscopic elec- 
tromagnetic stress tensor appropriate for ionized media. The general- 
ized Lorentz force includes the effects of polarization forces as well as 
those on the free charge and current densities. The resulting tensor 
is written in terms of the fields D, B, E, and H. Its expression for 
a fully ionized medium subject to an external electromagnetic field is 
discussed, as are the plasma conservation equations. An apparatus is 
suggested for its experimental discrimination. 



1 Introduction 

While the microscopic form of the electromagnetic stress tensor is well known 
PQ, its macroscopic form continues to be a matter of debate [21 [3], S] - Of inter- 
est are the effects of the material polarization fields P and M, particularly 
the magnetic contribution which has variously been given [3 |6] by Lorentz 
and Kelvin as F L k = /-to J x H + /i M ■ VH and by Korteweg and Helmholtz 
as F K h = J x B — H ■ HV/x/2. Following the arguments presented by 
Mansuripur [7J [8], we propose a form for the macroscopic electromagnetic 
stress tensor for ionized media expressed in terms of the fields D, B, E, and 
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H. As no assumption on the form of the constitutive relations is taken other 
their definition, we believe this expression has applicability beyond that of 
linear, homogeneous media. 

With restriction to a fully ionized medium P —>■ 0, we evaluate the non- 
linear relation between M and H explicitly. For vanishing charge density, 
the stress tensor may be expressed in terms of the pressure p, the magnetic 
field H, and the applied field E. The electromagnetic stress may then be 
combined with the hydrodynamic stress to yield the equation for the net 
conservation of momentum. With reduction to a stationary plasma carrying 
no free current in the absence of gravity, there remains a balance of pressure 
against the magnetic polarization force, suggesting a means for an experi- 
mental apparatus to determine whether its effect is present in a fully ionized 
medium. 

2 Field equations 

The macroscopic inhomogeneous Maxwell equations for the continuum (fluid) 
description of electromagnetic phenomena are 

V D = pi , V x H - d t D = J f , (1) 

where d t = d/dt and the sources appear opposite the fields. The correspond- 
ing homogeneous equations are 

V • B = , V x E + d t B = , (2) 

and we emphasize that the proper interpretation of these two equations [9|[T0] 
is not as dynamical equations but rather as geometric constraints [TT| [T2] 
yielding the potential formulation E = — V<P — dtA and B = V x A, and 
the gauge condition V ■ A + fi e dt^> = expresses the continuity of the 
potential [T3| associated with the continuity of the source V ■ J + d t p = 0. 
The constitutive relations are taken in their most general form, 

D-P = e E , B-M = fx H , (3) 

without assumption on their detailed expression. Note the units M = p, Msi 
used in this section and the next. 

In the following, the operator V is allowed to act upon the material 
polarizations P and M. Such allowance is consistent with the expressions 
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Vx P = Vx D and V-JVf/// = — V-i3" necessary to compute the extruded 
fields from the bar electret and bar magnet [2] . No requirement of linearity 
is imposed, as the tensor structure developed accommodates naturally any 
variation in the material polarizations. 

3 Generalized Lorentz force 

In its microscopic form, the Lorentz force density is written as 

i q 

where the sum is over the particle label q and the fields E and B are to be 
evaluated at the particle position. The most common macroscopic general- 
ization is 

F L = J2 n se s (E + V s xB)^pE + JxB, (5) 

s 

where the sum is over the species label s and the fields E and B are aver- 
aged over the infinitesimal volume (fluid) element containing sufficient par- 
ticle number for continuum thermodynamics to be applicable. An averaging 
process is inherent in taking the particle velocity v q over to the continuum 
velocity V s = ^2 ges v q /n s which here includes any diamagnetic gyration. 
The standard reduction gives the force density F L = V • T EB — e fi dtSEB in 
terms of the stress tensor T EB = e (EE - I E 2 /2) + fi l (BB - I B 2 /2) and 
the Poynting vector Seb —Ex B/hq = c 2 Keb for vacuum light speed c. In 
polarizable media [TJ [H] , one uses the Poynting vector Ser = E x H which 
accounts for the work done (energy deposited) to the free source density, 
suggesting that the most appropriate fields for the macroscopic description 
are E and H. 

We call the reader's attention to the excellent analysis by Mansuripur 
E], particularly the statement that "a new term had to be introduced 
to account for the force experienced by magnetic dipoles." This magnetic 
polarization force is familiar to anyone who has ever held two pieces of per- 
manent magnet close together, accounting for the attraction or repulsion. On 
energetic grounds [T5J ITS] , one expects an isolated dipole m to experience a 
force 

V(m • B) = m x (V x B) + (m • V)B (6) 
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which accounts for the interaction with both currents and fields, and similarly 
for the electric polarization. Mansuripur [5] has recently presented an 
analysis of the generalized Lorentz force for polarizable media in the absence 
of sources, 

i^mcdia = (P ■ V)E - (d t M) x e E + (M • V)H + (d t P) x p H , (7) 

and much of the following as regards the media contribution appears there 
first. To this expression we add the source contribution yielding the net 
Lorentz force = -F me dia + Efiee, where 

F free = p f E + J f x p H . (8) 

Why no J{ x M term? From Eqn. we see that such effect will be accounted 
for by the polarization term; furthermore, its additional presence would spoil 
the symmetry of the expressions below, hanging off the end like an unwanted 
tail. 

One proceeds by expressing -F me dia in terms of the material polarization 
tensors PE and MH using the identity V • (PE) = (V • P)E + (P • V)E 
to write 

fmedia = V • (PE) - (V ■ P)E - (d t M) x e E (9) 
+V ■ (MH) - (V • M)H + (d t P) x fi H , (10) 

then substituting d t P = V x H — J { — e d t E and d t M = —V x E — ^ d t H 
along with V • P = p { — e V • E and V ■ M = — yU V • H yields 

fmedia = V ■ (PE) + e [(V ■ E)E +(VxE)xE] (11) 
+V • (MH) + po [(V • H) H + (V x H) x H] (12) 
-e p d t (E x H) - Pi E - J f x p H . (13) 

Now adding the source term -Ff ree and expressing in dyads gives 

F L = V ■ (PE) + e [V ■ (EE) - VE 2 /2] (14) 
+V • (MH) +p [V ■ (HH) - VH 2 /2] (15) 
-e p d t (E x H) , (16) 

and one final substitution yields 

F L = V-(DE) + V-(BH)-V(e E 2 + p H 2 )/2-e p d t S ER (17) 
= V-T™-d t K m (18) 
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as the macroscopic electromagnetic stress tensor for source bearing polar- 
izable media in terms of the fields D, B, E, and H, noting that e^E 2 = 
(€qE) ■ E = (D — P) ■ E and similarly for fioH 2 . One may compare this 
expression with those by Abraham [17] and Minkowski [18J as well as the 
method by Louden et al [2J. This tensor form benefits from the verification 
of energy and momentum conservation by Mansuripur [7J, [S] in a variety of 
situations. 

Mansuripur considers also a media force which neglects the material po- 
larization tensors PE and MH in Eqn. (jUJ), 



n'edia = "(V • P)E - (d t M) x e E - (V • M)H + (d t P) x » H , (19) 



which Louden and Barnett [19] show yields an equivalent force given certain 

conditions on the material and may be expressed as V • T EH — d t K E n where 
T EH = Teh ~ ~ MH. Mansuripur [i8] gives the corresponding torque 
density expressions 



and we see the difference in the expressions for the torque about an observer 
r->0isPx£ + M x H, which may be related to the isolated dipole 
torques t p = p x E and t m = m x B. These torques serve to align the 
material polarizations along the fields in a fluid medium. 

4 Plasma magnetization 

Considering now one's favourite fully ionized plasma with species s G {e, i} 
and ion charge e» = ^e, we take P — > such that D = eoE, noting that the 
presence of neutrals and partials may contribute to an electric polarization 
P. With a total particle number of n = n e + rij, the electron density may 
be written n e = n — Pf/e, where n = riiZi is the number of electrons 
required for charge neutrality. Assuming equipartition of internal energy 
Wj~ = 2Ws = T s where T s <— fcsT s , the kinetic pressure is p = nT = 
J2 s n sT s and the mass density p m = nm = ^ s ra s m s , exhibiting the utility 
of distinguishing the expressions n and uq. For neutral, hydrogenic plasma, 
Zi — 1 and n = 2n . The free sources are defined by p f = J2 S n s^s and 
Jf = J2 s n s e s V s , and the fluid momentum by K { = p m Vf = J2 s n s m sV s - 






r x F mcdia + PxE + MxH, 



(20) 
(21) 
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From these relations, one may write V it€ = Vf±m eti Jf/ep m . The velocity V s 
here is understood not to include the diamagnetic gyration J dia = V x M$i 
which is associated with the diamagnetic momentum -RTdia- In terms of the 
gyro- vector oj s = -e s B s /m s , one writes K dia = p m V d ia = V x L dia for 
angular momentum Xdia = S s n s(^ s )> w h ere (h) = ^T s uj s /uj 2 . 

We now shift our notation so that the magnetic material polarization is in 
SI units M S i <— M / p and drop the subscript. The species dipole moment 
per unit volume is taken as M s = n s (m s ) = -n s (W^/B 2 )B s , where the 
field felt by a single particle of species s within the unit of volume is the 
net field less the particle's own contribution B s /p = H + M — (m s ) = 
H + M k + a s M s , where k ^ s and a s = (n s — l)/n s = 1 — l/n s is a unit-less 
factor. For a sufficiently dense plasma, a s — > 1 and B s — > B. The material 
polarization is given by the net dipole density [20], which for p = p/po and 
h = H/H may be written 

M = J2 Ms = Mme) +rn(mi) , (22) 

s 

s 

= -hM = -hp/ (H - M) , (24) 

and has the physical solution M/H = (1 — y/i — Ap/H 2 )/2 as the plasma 
is diamagnetic [21J. That simple expression represents a significant result, 
as we are aware of no other author investigating its utility. Ultimately, the 
proper treatment of magnetization requires the use of quantum theory, in 
particular as to account for spin [22] • 

From the form of the solution for M one can immediately read a limit 
on the ratio of kinetic to free magnetic pressure, (3h = 2p/H 2 < 1/2 for a 
dense plasma. In terms of the net field B, we have (3b = 2p/(H — M) 2 = 
Ph/0- ~ M/H) 2 < 2, and the ratio M/H is limited to 1/2. In the dilute 
fluid limit n Q — > 1/rf such that a = (n — l)/n = 1 — 1/2% — > 1/2, we 
find M s -> m s so that the limits (3 H -> M/H -> 1 when T { = T e . One 
must be careful to define the appropriate unit of volume for a dilute plasma, 
which we feel should be on the order of the cube of the ion gyro-radius 
Ti = Ut l ^j2Ti/mi. 

The magnetized plasma stress tensor for a fully ionized medium, T^h ~ 
Teh f° r -P — ► 0, may then be expressed as 

Tff = e (EE - lE 2 /2) + p [HH(l + y/l-Ap/H*) - IH 2 }/2 , (25) 
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where the electric field is the sum of that from plasma sources and any applied 
field E = E p + E app . In the neutral limit V • E p — > 0, the plasma electric 
field is dynamic, V x E p = —fiodt(l — M/H)H, thus specifying the tensor 
in terms of the pressure p and the fields H and E app . 



5 Conservation equations 

Most analyses of plasma [23J 12H ES] rely on the collisional Boltzmann equa- 
tion 

df df df (df\ fdf\ 

= C + S, (27) 

where C and S represent the collision and source terms whose moments C&, 
Sk are indexed by the natural numbers k G N . The species individually 
satisfy the particle continuity equation d t n s + V • n s V s = S s o where n s = 
[d s + (V • V s )]n s is the particle source rate and d s = d t + V s ■ V is the 
convective derivative, thus yielding the conservation of mass 

d tPm + V-K i = p m , (28) 

where p m = J2 S p s = J2 S "s m s is the mass source rate, and the conservation 
of charge arises naturally from Noether's theorem [26] applied to the gauge 
condition. 

The species contribution to the free momentum is K s = p s V s . For 
P — > the macroscopic Lorentz force reduces to 

F L = F free + p {M-V)H- Li e (d t M) x E , (29) 

and the species components may be extracted. With an influx S s \ ^ 
of mechanical momentum F s \ = dtK s i, the balance of source, acceleration, 
and force for each species K s = F s is written 

p s V s + p s {d t + V s ■ V)V S = F sl + F sG + F sL + F ss + F sk , (30) 

where F s q = — Vp s G accounts for gravitational acceleration in potential 
G, -F s l is the generalized Lorentz force on the s component of the fluid, 
the term F ss = — V ■ (p s + II S ) represents intraspecies collisions, and the 



7 



term F sk = Yls'^s F ss > ~ — p s v s k(V s — Vk) represents interspecies collisions. 
Using Fks = —F S }., their sum 

J2k s = f 1 + f g + f l - v-( P + n) , (31) 

s 

is the equation for the net conservation of momentum, which may also be 
written 

K { + V • (m e m t /e 2 p m ) J { J { = d t (K 1 - K EH ) - V • (p + II + IU G - Tff ) , 

(32) 

where U G = p m G is the gravitational potential energy density, making ex- 
plicit the appearance of the field momentum. The term in «7f«7f represents 
the free current's contribution to the convective force. 

Noting that the hydrodynamic and thermodynamic pressures may dif- 
fer by a factor of the bulk viscosity [27J, our identification of (the diver- 
gence of) the pressure and viscosity tensors as a collision term relies on 
an argument presented by Woods [23] that the mediator of the pressure 
force is the existence of scattering events and the observation that the Bra- 
ginskii viscosity [29] is written in terms of an energy density, a collision 
rate, and a shear tensor, n s = J2 a Vs^s f° r Vs ~ n sT s /vf s and W s = 
VF S + (W S ) T — (2/3)(V ■ V S )I. We would like to express the interspecies 
collision term F s k along similar lines in order to account for the momentum 
transfer in detailed form, postponing such development for a later time. The 
term for incoming force F\ = Y2 S Fsi — m s J v (dtf s )sd 3 v results from 
a source distribution (d t f s )s which is not isotropic in the particle velocity 
v s = V s + u s , thus yielding a net transfer of momentum. 

The kinetic energy density for each species is written 

8 S = Tr( Ps + n s + n s m s V s V s )/2 = 3 Ps /2 + n s m s V s 2 /2 , (33) 

where p s is the scalar pressure, whose evolution is described by the equation 

S s + V-q s = C s2 + S s2 + V S -F S} (34) 

in terms of the heat flux density q s , which may also be written [25] as 

d t S s + V-Q s = C s2 + S s2 + V s ■ (F s - F ss ) , (35) 

in terms of the energy flux density Q s = £ S V s + (p s + II S ) • V s + q s . Using 
J2s(C S 2 + V s ■ Fsk) = for elastic collisions, one may write the net energy 



8 



conservation equation 



+ = J2(c s2 + S s2 + V.-K.} , (36) 

s 

where the source moment Yl s gives the net influx of kinetic energy. Clo- 
sure is achieved by specifying the quantities q s , p s + II S , and u ss > in terms of 
the degrees of freedom n s , T s , and V s and the fields E and H . 

When evaluating the net kinetic energy of the plasma, one finds terms 
which depend on the current as well as the momentum. Using the notation 
Jf = Jf/ep m and defining S{ = 3p/2 + p m (V { 2 + m e miJf )/2 , one may write 
the rate of change of the net kinetic energy density as 

^ S a = S{ + V ■ 3(m e pi - mip e )J { /2 (38) 

s 

+ V ■ |m e mjp m Vf + (m e -mi)J{/2 -J f |j f , (39) 

where £f = d t £{ + V • £{Vf. Net energy conservation is best interpreted as 
the dynamical equation for the pressure, as it explicitly relates dtp to other 
quantities determined elsewhere. 



6 Ohm's law 

One final piece of the puzzle remains: determining the relationship between 
the electric field and the current, which in its simplest form reduces to the 
linear relationship Jf = oE. For a medium in which all species of charge car- 
riers may flow, the relationship is not so simple. The theory provides for two 
types of coupling, mass and charge, by which forces may act on particles, and 
while momentum is driven (along H) by forces acting on ions and electrons 
in the same direction, current is created by forces which drive them apart. 
(Perpendicular to H, life gets a bit more complicated.) Consequently, the 
generalized Ohm's law is given by (the sum of) the ion momentum conser- 
vation equation(s) minus the electron equation, which may also be thought 
of as an equation for positrons going backwards in time [15J. 

Many authors derive Ohm's law "by taking a particular linear combina- 
tion of the fluid equations" [23] where the ion equation is reduced by a factor 



9 



of m e /rrii, retaining only terms which survive the limit m e <C m^. We argue 
that such procedure does not tell the whole story, effectively neglecting the 
remainder of the ion contribution (1 — m e /rrii), nor does it respect the un- 
written factor of units inherent in physical equations (as opposed to purely 
mathematical ones). Before adding or subtracting two physical quantities, 
they must be expressed in the same units, including prefix. The extrane- 
ous mass factors applied when taking m e Fi — m,iF e may just as well be 
replaced by "milli" and "mega" as prefix on the originally common units, 
so that the difference has numerically subtracted mega-Newtons from milli- 
Newtons. Just as the free momentum equation relies on the net force balance 
Ki + K e = Fi + F e , so does the free current equation rely on the net force 
difference Ki — K e = Fi—F e , which through substitution yields the equation 
for Jf. From previous derivation [30J, we note that the extraneous factors 
institute a simplification of the equation through the artificial cancellation of 
certain terms, particularly in the convective derivative, thus revealing their 
likely motivation. 

7 Discriminatory apparatus 

After so much development, we restrict ourselves to a single evaluation. Con- 
sider an annular plasma chamber in (R, (p, Z) coordinates surrounding a con- 
ductive wire along the Z-axis which is insulated from the plasma, as in Fig. [TJ 
Supposing the plasma to be stationary and carrying no free current, the net 
force balance equation reduces to 



in the absence of gravitational acceleration g, where we have used d t = 
and V x H = in the plasma region. Without the media contribution to 
the Lorentz force, one has simply Vp = thus a constant pressure in the 
chamber, and the current in the wire should have no effect on the plasma, 
thus offering a chance for an experimental determination of its existence. A 
thorough analysis would need to account for the momentum transfer by the 
electric field during the period when the current is changing. 

Sufficiently far from the ends of the wire, its field is given by Ampere's 
law as H^ irc = (f)H /R for H = Iz/2tt, where Iz is the total current, and 



Vp = // (M • V)H 



IM>(M/H)(H-V)H , 
fx (M/H)VH 2 /2 , 



(40) 
(41) 
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its magnitude has the gradient Vif 2 /2 = —RH^/R 3 . Using p = p/fio, the 
net force balance becomes the scalar equation 

dRP = -\(l-^jl-4v/Hl)d R Hl, (42) 
= l^-^^RVBl^H^R- 3 , (43) 

which one may solve for the pressure profile given values for the current Iz 
and the pressure at the outer chamber wall p c . 

For a plasma chamber with R c = 100 mm outer radius and inner radius 
of 5 mm surrounding a current of 50 A along Z at R = 0, we calculate the H 
field profile in Fig. [2](a). Assuming a common temperature of 100 eV and a 
net edge density of 2 x 10 13 m~ 3 , the pressure profile in Fig.[2](b) is normalized 
by the pressure at the chamber wall p c = 0.32 mPa. In the presence of 
-Fmedia, the magnetic polarization force is balanced by a pressure gradient 
such that the pressure at the centre is reduced by a factor ~50 relative 
to its value at the edge. Without that force, there should be no pressure 
gradient. The location of maximum (3h is at the outer wall R c , and for these 
parameters it remains below its limit of 0.5 as shown in Fig. [2]^c). The upshot 
is that experimentalists should be able to construct an apparatus which can 
measure the effect predicted by the theory, requiring verification that fluid, 
diamagnetic material should be repelled from a stronger field region. 

8 Conclusions and outlook 

In summary, we have found that the Lorentz force on the free source densi- 
ties may be incorporated self-consistently with the media force proposed by 
Mansuripur [7J [8] and investigated by Louden et al [21 [19] to yield the macro- 
scopic electromagnetic stress tensor in terms of the fields D, B, E, and H 
without requirements on the detailed form of the constitutive relations. The 
nonlinear model for plasma magnetization is given a solution which displays 
a limit on the ratio of kinetic to magnetic pressure. The macroscopic con- 
servation equations for a fully ionized medium are expressed in terms of the 
free momentum and current, making explicit the contribution of each to the 
kinetic energy density. 

With restriction to a stationary plasma without free current and neglect- 
ing gravity, conservation of momentum requires a pressure gradient develop 



11 



to oppose the magnetic polarization force. An evaluation for parameters 
which may be accessible to experimentalists suggests a measurable difference 
between the central and edge pressures for a chamber surrounding a steady 
current, indicating the presence of the media contribution to the generalized 
Lorentz force. 

The complete form of the macroscopic Lorentz force bears directly on 
the physics of plasmas. Treatments by the fluid description need to respect 
the dielectric and diamagnetic properties of the medium as encoded in the 
constitutive relations. The model presented here is not the only one available, 
and there is something to be said [31] for using the expression [V x (JEf + 
M)] x (H + M) + V(iW • H) as the magnetic contribution to the Lorentz 
force. Variations in a theoretical development are no bad thing, as they allow 
for the possibility to discriminate between competing models. We consider 
the full effect of plasma magnetization to be an open question, requiring a 
conspiracy between theorists and experimentalists for its determination. 
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Figure 1: Plasma chamber surrounding a current carrying wire. A pressure 
gradient should develop to oppose the magnetic polarization force. 
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Figure 2: Solution for a plasma chamber with outer radius of R c = 100 mm 
and inner radius of 5 mm surrounding a wire at mm carrying 50 A current 
upward, (a) Magnetic field profile, (b) Normalized pressure profile for p c = 
0.32 mPa. (c) The (3h profile remains below the limit of 0.5 at the edge of 
the chamber. 
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